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Abstract 



In this paper we briefly review a model that describes the diffusion- 
controlled aggregation exhibited by particles as they are deposited on a sur- 
face. This model allows to understand many experiments of thin film deposi- 
tion. In the first part, we describe the model, which incorporates deposition, 
particle and cluster diffusion, and aggregation. In a second part, we study 
the dynamical evolution of the model. Finally, we analyze the effects of small 
cluster mobility, and we show that the introduction of cluster diffusion dra- 
matically affects the dynamics of film growth. Some of these effects can be 
tested experimentally. 
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I. INTRODUCTION 



Understanding the processes underlying the growth of thin films has led to widespread 
interest [jl|, both from theoretical and experimental points of view ||^12|. Recently, 
improvements in experimental techniques — such as scanning tunneling microscopy — permit 
the investigation of atomic details of the embryonic "sub-monolayer" stages of film growth. 
Recent experimental works |[TD|-|T^ have shown that fractal structures can be formed during 
the first stages of the growth. Then, it is interesting to try to analyze the growth in the 
framework of fractal models. One might consider the use of the percolation model [0,0 to 
describe certain experiments of surface deposition [|^. However, percolation assumes that 
particles do not diffuse after being deposited, when in fact not only diffusion but also aggre- 
gation of the diffusing particles takes place. There exist models of diffusing particles that 
aggregate, but such "cluster-cluster aggregation" (CCA) models |jT6| do not allow the con- 
tinual injection of new particles via deposition. Here we develop a model that incorporates 
the three physical ingredients of thin film growth: deposition, diffusion and aggregation 
(DDA). Similar models that neglect the shape of the islands or the possibility of cluster 
diffusion were studied independently 0J^. We have shown in detail elsewhere p how the 
DDA model generates a wide variety of fractal structures characteristic of different models 
such as percolation, diffusion limited aggregation (DLA) ]T7| , pr8| or CCA. Here we focus on 
the dynamics of film growth and the importance of (small) island diffusion. 



II. MODEL DESCRIPTION AND JUSTIFICATION 

The DDA model is defined as follows (Fig. H): 

(1) Deposition. Particles are deposited at randomly-chosen positions of the surface at a 
flux F per lattice site per unit time. 

(2) Diffusion. All particles and clusters (sets of connected particles) are chosen at random 
and attempted to move North, East, South or West by one lattice constant per unit time. 
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The probability that they actually move is proportional to their mobility, which we assume to 
be given by Dg = Dis^^' . Here s is the number of particles in the cluster, Di is the diffusion 
coefficient for a monomer (s = f ), and the parameter 7 characterizes the dependence of Dg 
on cluster size. 

(3) Aggregation. If two particles come to occupy neighboring sites, they (and therefore 
the clusters to which they belong) stick irreversibly. 

We call particles the isolated atoms (or monomers) that are deposited on the surface, 
clusters any set of connected particles (including the monomers) and islands the clusters 
containing more than one particle. Physically, two competing mechanisms are introduced 
in the model, each one with its own time scale: deposition and diffusion. It is useful to 
introduce the normalized flux defined as the number of particles deposited per unit site per 
diffusion time r, where r is the mean time needed by a monomer to jump by a lattice site. 
The monomer diffusion coefficient is then given by Di = l/(4r), and the normalized flux 
by = Ft. Then, from experimental values of F and Di it is possible to calculate (p and 
the morphologies predicted by our model. The program actually calculates a probability for 
dropping a particle: Pdrop = 4>L'^/ (0-^^ + N^i) where L is the system size and is the total 
number of clusters present in the system. A random number p is chosen and compared to 
Pdrop- li p < Pdrop, a particle is added at a random position on the lattice. If p > pdrop, a 
cluster or a particle is chosen at random and attempted to move. In both cases, the time is 
increased by r/ (0-L^ + N^). 

It should be stressed that this is only a "zeroth-order" model which has the ambition to 
give a feeling on the relative influence of deposition and diffusion on the growth properties of 
films. Details specific to certain experimental systems, such as the existence of the Schwoebel 
barrier, the precise dependence of cluster diffusion on size, etc. are not carefully taken into 
account since we want to keep the DDA model as general as possible (see Ref. 0). 
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III. DYNAMICAL EVOLUTION 



We present here snapshots of the system at different times to show the formation of the 
islands. These images were obtained for = 1.2 10^^ and 7 = cxo (i.e., only monomers are 
allowed to move). A detailed presentation of the dynamical evolution of the model has been 
given elsewhere 0. Actually, it is more interesting to use the surface coverage 9 instead of 
the time as the evolution parameter. 9 is defined as the ratio of the number of occupied sites 
to the total number of sites on the surface, and for the times studied here we have 9 ~ Ft. 

At very short times {9 < 0.001), mainly monomers (isolated particles) are found on 
the substrate, since they did not yet met another one to form a cluster. Later, small 
clusters are homogeneously grown on the surface (Fig. |^a), and the island density (i.e. the 
number of islands per lattice site) starts to grow. These small clusters can be considered 
as the "nucleation centers" for the growth. As time increases, large clusters grow on these 
nucleation centers, by addition of single particles (Fig. Qd). These clusters are very similar 
to those obtained experimentally (Figs. |l|a-b). We have found that their fractal dimension 
is 1.65, very close to the fractal dimension of the DLA clusters. This (DLA-like) growth 
mechanism goes on until the linear dimension of the clusters becomes comparable to the 
separation between them (Fig. ^). Then, many particles start to fall inside the clusters 
and their fractal dimension rapidly increases. Eventually a cluster of a size comparable to 
the system size is built and the system spans (Fig. ^). At that time, the effective fractal 
dimension of the spanning cluster reaches a value close to 1.9. 

It is interesting to note that the growth of the film for 7 = 00 is in some loose sense "self- 
regulated" since islands tend to avoid each other, by growing preferentially in the directions 
where no island is present. The reason is that monomers are more abundant in those regions 
since no island captures them. Then, the probability for a given point of an island to grow 
is higher if it is far from other islands. This effect can be quantified in the following way. 
We first grow a system with a given coverage and we freeze the growth. Then, we deposit a 
monomer and we move it exactly as during the growth, but each time the monomer reaches 
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an island, we record the point where it touched the island and remove it. By counting the 
number of monomers that hit each point, we can measure the growth probability within each 
island. Fig. ^ shows a global view of these probabilities for a large system at a coverage 
of .04. The white points correspond to those point which grow slower and the bright blue 
points show the regions of high growth probability. Fig. ^ shows in detail the growth 
probabilities for one of the islands. The fact that each island has its own "capture zone" 
and grows inside it has some interesting consequences on the size distribution of the islands. 
These will be discussed below. 



IV. ISLAND DIFFUSION 

In the last paragraph we have studied the behavior of the model in the case when only 
monomers are allowed to move. In this section, we want to address some consequences of 
(small) cluster diffusion. It has long been recognized that cluster diffusion can influence the 
growth of the films 0,|19|, even if at that time experimental proofs were lacking. Recently, 



experiments [^] and molecular dynamics studies have shown that small clusters can move on 
the surfaces without breaking |2l|. Generally, however, the experimental results are analyzed 



in the framework of the different models containing only monomer diffusion. These models 
can include reversible aggregation meaning that two monomers that aggregate can 



detach after a certain time. This can lead to wrong interpretations if small cluster mobility is 
indeed present in the experiments [|23|. Then, it is important to understand the consequences 
of this mobility to be able to identify them in the experiments. 

We have shown previously that the introduction of cluster mobility considerably 
changes the growth dynamics of the film. Specifically, we showed that: 

(z) The mean cluster size increases exponentially as a function of the coverage. A power- 
law dependence is generally found when only monomer diffusion is allowed. This effect is 
intriguing and may be due to large cluster diffusion 0. 

(m) The maximum island density depends on the incident flux according to Nmax ~ -^^ 



5 



with ( = 0.42, as opposed to ( = 0.33 found with monomer diffusion. It is also known 
that when only dimer diffusion is introduced, C = 0.4 ||^. The problem is that generally 
the accuracy of the experimental determinations of Nmax is not high enough to distinguish 
between the different exponents |]12| . 



{in) More interesting experimentally: the evolution of the island density as a function 
of coverage for low normalized fluxes (typically less than 10^^). It is known [|] that the 
maximum of the island density is reached only at roughly a coverage of 0.2 when only 
monomers can move. On the contrary, if clusters can move, this maximum is reached for 
smaller values of the coverage. 

We add here a test that is easy to perform experimentally, and which has actually already 
been used to distinguish different growth mechanisms [Il2| , p3| , p4| : the rescaling of island size 
distributions @J2^. The idea is to rescale the island size distributions into a universal 
distribution that depends on the ingredients of the model (i.e. the detailed mechanisms of 
Deposition, Diffusion and Aggregation) but not on the values of the flux or the coverage. 
Very recently, Mulheran and Blackman have given some interesting insights of why 
such a universal function should exist, at least for the case of heterogeneous growth. Their 
argument can be summarized as follows. At the beginning of the growth, nucleation centers 
form (Fig. |^a). Then, each center grows by catching the monomers falling inside its "capture 
zone", roughly identified with its Voronoi polyhedron. Therefore its size is, at any time, 
proportional to the surface of its Voronoi polyhedron, which does not change with time (if 
one neglects nucleation of new islands in the case of homogeneous nucleation). The result is 
that at any coverage the size distribution of the islands reproduces that of the Voronoi cells, 
which explains the rescaling for different coverages. We show in Figure ^ that the scaling 
of the size distributions also occurs when clusters are allowed to move, but the universal 
function is different from that found when only monomers move. Three different diffusion 
hypothesis have been made: only monomers move, monomers and dimers diffuse, all clusters 
up to size 100 do move (7 = 1 has been taken in the two last cases). We see that the universal 
function becomes narrower and has a higher maximum when larger and larger clusters are 



allowed to diffuse. 



V. DISCUSSION AND CONCLUSION 

In summary, we have proposed a model that describes the diffusion-controlled aggrega- 
tion exhibited by particles as they are deposited on a surface. The model, which incorporates 
deposition, particle and cluster diffusion, and aggregation closely reproduces some exper- 
imental images (compare Figs [l|a,b and Figs We find that the model permits one to 
distinguish the effects of deposition, diffusion and aggregation, and that tuning the relative 
strength of, e.g., deposition and diffusion, generates a rich range of morphologies — including 
diffusion limited aggregation, cluster-cluster aggregation, and percolation. The length and 
time scales characterizing these morphologies depend on such experimentally-controllable 
parameters as deposition flux and diffusion constant, raising the possibility that the model 
may prove useful in future studies seeking the controlled design of nanostructure morpholo- 
gies. We can argue that the DDA model is suited for (^) MBE by taking large 7 values, (m) 
Other deposition experiments where edge diffusion is absent, for example aggregate deposi- 
tion |jl2[|. {in) All the dynamic properties related to sma// cluster diffusion on surfaces such 
as those presented in Sec. VI. 
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FIGURES 

FIG. 1. Experimental images of the first stages of growth for films prepared by (a) Atomic 
deposition of Ag on Pd [0| and (b) Deposition of Sb compact aggregates containing 2300 atoms 



on graphite |12|. 



FIG. 2. Schematic representation of the basic processes considered in this model: (a) deposition, 
(b) and (d) particle diffusion, (e) island diffusion and (c) aggregation, (b) corresponds to nucleation 
(i.e. a new island (c) is created) while (d) corresponds to growth of an already existing island (see 
the text for details). 

FIG. 3. Morphologies obtained for a normalized flux (j) = 1-2 10~^, a system size L = 500 and 
7 = oo. The images show a portion 290 x 210 of the lattice. Four different coverages (corresponding 
to four different times of deposition) are shown: (a) coverage 0.02 (b) coverage 0.1 (c) coverage 
0.25 (d) coverage 0.4. Colors indicate the order in which particles were deposited on the surface: 
the first particles deposited are colored white, then yellow, orange, red, blue and finally green. 

FIG. 4. Growth probabilities for different points of the islands at a coverage 0.04 (7 = 00). 
(a) global view of the surface (b) detailed view of an island. Colors give the probability that a 
given site of the island has to catch a monomer. The probability increases in the order: white, 
yellow, orange, red, blue, green and finally big bright blue dots. 

FIG. 5. Rescaled island size distributions. The island distributions are transformed to p(s) 
which represents the probability that a randomly chosen cluster belongs to an island containing s 
clusters The island size s is scaled by the mean island size Sm- Sets (1) to (3) correspond to 
simulations where only monomers are allowed to move (7 = 00). For sets (4) to (6), 7 = 1 but 
only dimers can move, while for sets (7) and (8) clusters containing up to 100 monomers can move. 
The detailed parameters used for each set are as follows: (1) 4> = 10~®, coverage 0.3 ;(2) = 10~^, 
coverage 0.15 ; (3) i;^ = 10~^, coverage 0.1 ; (4) = 10~^, coverage 0.05 ; (5) = 10~^, coverage 
0.1 ; (6) = 10^^, coverage 0.3 ; (7) </> = 10"^, coverage 0.15 ; (8) (/> = 10"*^, coverage 0.05 ; (9) 
(j) = 10~^, coverage 0.1. The lines represent averages of the distributions for each case. 
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